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Abstract 

If the atmospheric neutrino oscillation amplitude, sin^ 2^atm is not maximal, there is 
^ a two fold ambiguity in the neutrino parameter space: sin^ ^atm > 0.5 or sin^ 6'atm < 0.5. 

O ' In this article, we study the impact of this degeneracy, the so-called octant degeneracy, 

on the T2KK experiment, which is a proposed extension of the T2K (Tokai-to-Kaimoka) 
neutrino oscillation experiment with an additional water Cerenkov detector placed in 
^1 Korea. We find that the degeneracy between sin^ ^atm = 0.40 and 0.60 can be resolved 

' at the 3(T level for sin'^ 20rct > 0.12 (0.08) for the optimal combination of a 3.0° off-axis 

beam (OAB) at SK (L = 295km) and a 0.5° OAB at L = 1000km with a far detector of 
^ ' lOOkton volume, after 5 years of exposure with 1.0 (5.0) x lO^^POT/year, if the hierarchy 

I is normal. We also study the influence of the octant degeneracy on the capability of T2KK 

experiment to determine the mass hierarchy and the leptonic CP phase. The capability 
of rejecting the wrong mass hierarchy grows with increasing sin^ ^atm when the hierarchy 
is normal, whereas it is rather insensitive to sin^ ^atm for the inverted hierarchy. We also 
find that the la allowed region of the CP phase is not affected significantly even when 
the octant degeneracy is not resolved. All our results are obtained for the 22.5 kton 
Super-Kamiokande as a near detector and without an anti-neutrino beam. 



1 Introduction 

A decade ago, it was difficult to believe that neutrinos have mass and the lepton flavor mixing 
matrix, the Maki-Nakagawa-Sakata (MNS) matrix [1], has two large mixing angles [2]- [7]. 
Within the three neutrino framework, 2 mass-squared differences, 3 mixing angles, and 1 CP 
phase can be resolved by neutrino oscillation experiments. So far, the magnitude of the larger 
mass-squared difference, the magnitude and the sign of the smaller one, two of the three mixing 
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angles, and the upper bound of the third mixing angle have been known. The sign of the larger 
mass-squared difference (the mass hierarchy pattern), the magnitude of the third mixing angle 
(^rct), and the leptonic CP phase ((5mns) arc yet to be measured. 

In the previous papers [8, 9], we studied in detail the physics impacts of the idea [10] of 
placing a far detector in Korea along the T2K neutrino beam line. For concreteness we examined 
the effects of placing a lOOkton water Cerenkov detector in Korea, about L = 1000km away 
from J-PARC (Japan Proton Accelerator Research Complex) [11], during the T2K (Tokai-to- 
Kamioka) experiment period [12], which plans to accumulate 5 x 10^^ POT (protons on target) 
in 5 years. We find that this experiment with two detectors for one beam, which may be 
called the T2KK experiment [13], can determine the mass hierarchy pattern, by comparing 
the i^fj, Ve transition probability measured at Super-Kamiokande (SK) and that at a far 
detector in Korea. Moreover, both the sine and cosine of the CP phase can be measured 
from the energy dependence of the — ^ oscillation probability, which can be measured 
by selecting the quasi-elastic charged current events. By studying these physics merits of the 
T2KK experiment semi-quantitatively, we find an optimal combination of a 3° off-axis beam 
(OAB) at SK and a 0.5° OAB in the east coast of Korea at L = 1000km, for which the 
mass hierarchy and the CP phase (^mns) can be detemined without invoking an anti- neutrino 
phase [8, 9], when the mixing angle ^^ct is not too small. In the related study [14], a grander 
prospect of the T2KK idea has been explored, where two identical huge detectors of several 
lOOkton volume is placed in Kamioka and in Korea, and the future upgrade of the J-PARC 
beam intensity has also been considered. The idea of placing two detectors along one neutrino 
beam has also been explored for the Fermi Lab. neutrino beam [15]. 

In this report, we focus on the yet another degeneracy in the neutrino parameter space, 
which shows up when the amplitude of the atmospheric neutrino oscillation, sin^2^ATM, is not 
maximal. Hereafter, we call this degeneracy between sin^ ^atm > 0.5 and sIu^^atm < 0.5, or 
between "90° — ^atm" and "^atm", as the octant degeneracy [16]. 

Since the best fit value of the mixing angle ^atm is 45° [2, 3, 4], we set sin^ 2^atm = 1 in 
all our precious studies [8, 9], and hence we did not pay attention on physics impacts of the 
octant degeneracy. However, we are concerned that the octant degeneracy affect the capability 
of the T2KK experiment for the mass hierarchy determination and the CP phase measurement, 
because the leading term of the z/^ oscillation probability is proportional to sin^ ^atm, not 
sin^ 2^ATM- If the value of sin^ 2^atm is 0.99, which is 1% smaller than the maximal mixing, the 
value of sin ^atm is sin ^atm — 0.45 or 0.55, which differ by 20%. 

For sin2 2^ATM = 0.92, which is stiU allowed at the 90% CL [2, 3, 4], sin^ ^atm = 0.64 or 
0.36, which differ by almost a factor of two. Therefore, we also examine impacts of varying 
sin^ ^^ATM on the mass hierarchy determination and the CP phase measurement by T2KK. In 
our semi-quantitatively analysis, we follow the strategy of ref. [8, 9] where we adopt SK as a 
near side detector and postulate a 100 kton water Cerenkov detector at L = 1000km, and the 
J-PARC neutrino beam orientation is adjusted to 3.0° at SK and 0.5° at the Korean detector 
site. 

This article is organized as follows. In section 2, we fix our notation and show how the octant 



degeneracy affects the oscillation probabilities. In section 3, we review our analysis method and 
present an explicit form of the Ax^ function which we use to measure the capability of the T2KK 
experiment semi-quantitativcly. In section 4, we show the results of our numerical calculation 
on the resolution of the octant degeneracy. In section 5, we examine the capability of the T2KK 
experiment for the mass hierarchy determination, in the presence of the octant degeneracy. We 
also show the effect of the octant degeneracy on the CP phase measurement in section 6. In 
the last section, we summarize our results and give discussions. 

2 Oscillation formular and experimental bound 

When a neutrino of flavor a is created at the neutrino source with energy E, it is a mixture of 
the mass eigenstates, u- 
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ka) = Z!^mki)> (" = e,/x,r) (1) 

i=l 

where Uai is the element of the Maki-Nakagawa-Sakata (MNS) matrix [1]. Without loosing 
generality, we can take Ue2 and to be real and non-negative and allow Ugs to have a 
complex phase 6mns [17, 18]. 

After travehng the distance L in the vacuum, a neutrino flavor eigenstate is found with 
the probability 
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Equation(2) shows that neutrino flavor oscillation is governed by the two mass-squared differ- 
ences and the lepton number conserving combinations of the MNS matrix elements. 

We take IA13I > IA12I without loosing generality. Under this parameterization, atmospheric 
neutrino observation [2] and the accelerator based long baseline (LBL) experiments, K2K [3] 
and MINOS [4], which measure the i/^ survival probability, are sensitive to the magnitude of 
the larger mass-squared difference and U^^: 

1.5 X 10-=^eV^ < \ml - m?| < 3.4 x 10-=^eV^ , (4a) 
sin^ 2^^™ = 4^723 (1 - U%) > 0.92 , (4b) 



each at the 90% confidence level. Hereafter, we use sin^ ^atm instead of the for brevity. 

The reactor experiments, which observe the survival probability of the Pg at L ~ 1km from 
a reactor, arc sensitive to the value of the larger mass-squared difference and the absolute value 
of Ue3- The CHOOZ experiment [5] reported no reduction of the Pg Aux and find 

sin2 2^i,cT = 4 \U,sf (1 - \Ue3f) < (0.20, 0.16, 0.14) 

(2.0, 2.5, 3.0) X 10- (5) 



for 
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at the 90% confidence level. In the following, we denote |C/e3p as sin^^ncT- 

The solar neutrino observations [6], and the KamLAND experiment [7], which measure 
the survival probability of and Pg, respectively, are sensitive to the smaller mass-squared 
difference and the value of Ue2- The present constraints can be expressed as 

ml-ml = (8.0 ± 0.3) x lO^^eV^ , (6a) 
sin^^soL = 0.30 ±0.03. (6b) 

The sign of ml — ml is determined by the matter effect in the sun [19, 20]. In these experiments, 
the order of A12 is roughly 1 and the terms with A13 oscillate quickly within the experimental 
resolution of L/E. After averaging out the contribution from A13, and neglecting terms of 
order sin^^^cT, we obtain the relation; 

sin^ 2^so. = ^U'^iU!, = AUUl - Ul^ - \ul^\) . (7) 

These simple identification, eqs.(4b), (5), and (7), are found to give a reasonably good descrip- 
tion of the present data in dedicated studies [21] of the experimental constraints in the three 
neutrino model. In this paper, we parameterize the CP phase as [18] 

(^MNs = - arg Uez . (8) 

The other elements of the MNS matrix can be obtained form the unitary conditions [17]. 
This convention allows us to express the MNS matrix directly in terms of the three observed 
amplitudes. 

The probability of the neutrino oscillation, eq.(2), is modified by the matter effect [19, 20], 
because only and Pg feel the potential by the extra charged current interactions with the 
electron inside the matter. This extra potential for is written as 

a = 2V2G,E^n^ ^ 7.56 x 10-^[eV^] (^) , (9) 



where is the Fermi coupling constant, E^, is the neutrino energy, 77e is the electron number 
density, and p is the matter density. The extra potential for Pg has the opposite sign. Because 
the matter effect is small at low energies and also because the phase factor A12 is small near 



the first oscillation maximum, A13 ~ vr, we find that an approximation of keeping terms linear 
in the matter effect and A12 is useful for analyzing the LBL experiments at sub GeV to a few 
GeV region [8, 9, 22, 23]: 



P,^^,^ = 2(1 + q) sin^ ^HCT (1 + A'') sin^ + S^) , (10a) 

P.^^.^ = 1 - (1 - g^) (1 + A^) sin^ + B^^ , (10b) 

where 

sin^ 6'atm = , (11) 

and A" and 3°" are the correction terms to the amplitude and the oscillation phase, respectively. 
For q; = e, we find 



A = -r — - COS 26'rct ^-^ \ — — sm b^^s , (12a) 

Ai^E 2 sm 6Irct V 1 + ? 



-B = -tt; cos26'rct + 77^-^^ Wt— cos^mns - sm 6'sol ■ (12b) 

4E 2 V 2sm6'RCT V 1 + ? / 



The octant degeneracy between "^atm" and "90° — ^^atm" corresponds to the degeneracy in the 
sign of q. When q denotes the true value for the octant degeneracy, —q is its fake value. 

Using typical numbers of the parameters from the atmospheric neutrino observation and 
LBL experiments, cq.(4), and those from the solar neutrino observation and the KamLAND 
experiment, eq.(6), the —>■ Vf. transition probability can be expressed as 



P.^_^ ~ 0.05 (1 + q) {^^^ (1 + A^) sin^ (- 
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around the oscillation maximum, | A13I ~ tt. Since the amplitude is proportional to sin^ ^^atm = 
(1 + q')/2, we expect that the octant degeneracy can be solved by measuring the z/^ — > z/g 
transition probability, i/the value of the sin^ 2^^pcT is known precisely. Because the first term 
of A^ changes sign according to the mass hierarchy pattern, A13 ~ vr for the normal and 
^13 ~ for the inverted, the amplitude of the transition probability is sensitive to the mass 
hierarchy pattern. The difference between the two hierarchy cases grows with the baseline length 
when L/ E is fixed at around the oscillation maximum [8, 9]. If there is only one detector at 



L ~ O(100)km, the small difference from the matter effect can be absorbed by the sign of q in 
the leading term of eq.(13a). 

The g-dependence in and in eqs.(13b) and (13c) may seem to affect the measurement 
of the leptonic CP phase. We find, however, that the g-dependence of the coefficient of the CP 
phase in the — > transition probability is not strong, because 
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which is independent of the octant degeneracy. Our numerical studies presented below confirms 
the validity of the above approximations. 

Around the first dip of the survival probability | A13I ~ tt, we find 
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Although the shift in the amphtude A^ and that in the phase B^^ are both proportional to q, 
their magnitudes are found to be less than 0.7% and 0.5%, respectively, for \q\ < 0.28, eq.(4b). 
Our numerical results confirm that the measurement of the i/^j, survival probability does 

not contribute significantly to the resolution of the octant degeneracy. On the other hand the 
smallness of the deviation from the leading contribution allows us to constrain |m| — | and 
sin^ 2^ATM accurately by measuring the — > i/^ survival probability. 



3 Analysis method 

In this section, we explain how we treat signals and backgrounds in our numerical analysis, 
and introduce a function which measures the capability of the T2KK experiment semi- 
quant it at ively. We consider a water Cerenkov detector at Korea in this study, because it allows 
us to distinguish clearly the events from jx^ events. The fiducial volume of the detector 
placed at Korea is assumed 100 kton, which is roughly 5 times larger than that of SK, 22.5 
kton, in order to compensate for the longer base-line length. We use only the CCQE events 
in our analysis, because they allow us to reconstruct the neutrino energy event by event [3]. 
Since the Fermi-motion of the target nucleon would dominate the uncertainty of the neutrino 
energy reconstruction, which is about 80 MeV [3], we take the width of the energy bin as 



6E^ = 200 MeV for > 400 MeV. The signals in the i-th energy bin, El = (200MeV x i) < 
E^, < EI + SE^, are then calculated as 

=MNa , <^u,{E) P.^^.^E) a^'^iE) dE , (16) 

JEl 

where P^^^i,^ is the neutrino oscillation probability including the matter effect, M is the de- 
tector mass, A^^ = 6.017 x 10^^ is the Avogadro constant, is the z/^ flux from J-PARC 
[24], and is the CCQE cross section per nuclcon in water [3]. For simplicity, the detection 
efficiencies of both detectors for both and Ve CCQE events are set at 100%. 

We consider the following background events for the signal of e-like events {a = e) and /x-like 
events (a = /i), 

Ar^^« = N'M) + Km + Ni{9,) , (a = e , /.) , (17) 

respectively. The three terms correspond to the contribution from the secondary neutrino flux 
of the primary beam, which are calculated as in eq.(16) where (^y^iyE) is replaced by ^^^{E) 
for 1^13 — ,T?e All the primary as well as secondary fluxes used in our analysis are obtained 
from the web-site [24] . After summing up these background events, the e-like and /i-like events 
for the i-th bin are obtained as 

K = Ki'^,) + K'''' , (« = e,/x), (18) 

respectively. 

Our interest is the potential of the T2KK experiment for solving the octant degeneracy and 
its influence on the resolution of the other degeneracies. In order to quantify its capability, we 
introduce a function, 

= XSK + xL + Xsys + Xpara > (19) 

which measures the sensitivity of the experiment on the model parameters. The first two terms, 
XsK Xkd measure the parameter dependence of the fit to the SK and the Korean detector 
data, respectively, 

where Nj^^^ is the calculated number of events in the i-th bin, and its square root gives the 
statistical error. Here the summation is over all bins from 0.4 GcV to 5.0 GcV for A^"^, 0.4 GeV 
to 1.2 GcV for Ne at SK, and 0.4 GeV to 2.8GcV for A^^ at Korea. In this energy region, we 
can include the second peak contribution in our analysis at Korea. We include the contribution 
of the yU-like events in order to constrain the absolute value of A13 strongly in this analysis, 
because a small error of A13 dilutes the phase shift B'^ [8, 9, 22]. 



A^f* is calculated by allowing the model parameters to vary freely and by including the 
systematic errors. We take into account four types of the systematic errors in this analysis. 
The first systematic error is for the uncertainty in the matter density, for which we allow 3% 
overall uncertainty along the baseline, independently for T2K (/p^) and the Tokai-to-Korea 
experiment {f^^)'- 

pT = f?l^r' P = SK, Kr). (21) 

The second ones are for the overall normalization of each neutrino flux, for which we assume 
3% errors, 

= 1±0.03, (22) 

for {vfj = i/g, Pg) i^/i)) which are taken common for T2K and the Tokai-to-Korea experiment. 
The third ones are for the CCQE cross sections, 

K*)"-/?^-?')"'". (23) 

where a denotes i = e — /i and £ = e — p,. Because i/g and i/^ CCQE cross sections are expected 
to be very similar theoretically, we assign a common overall error of 3% for i/g and and an 
independent 3% error for z/g and CCQE cross sections. The last one is the uncertainty 
of the fiducial volume, for which we assign 3% error independently for T2K (/y^) and the 
Tokai-to-Korea experiment (fy^). A^^'^* is then calculated as 

and accordingly, Xsys ^^^^ terms; 

To put them shortly, we account for 4 types of uncertainties which are all assigned 3% errors: 
the effective matter density along each base line, the normalization of each neutrino flux, the 
CCQE cross sections for z/; and and for the fiducial volume of SK, and that of the Korean 
detector. In total, our Ax^ function depends on 16 parameters, the 6 model parameters and 
the 10 normalization factors. 

Finally, Xpaxa accounts for external constraints on the model parameters: 



^P^™ [ 0.6 X 10-5 j ^ V 0-07 

^ |^ sin^2eT-^in^2fe"^ j ' _ (26) 




- mlf - 8.2 X 10-W\ ^ /sin2 2^^*^ - 0.83 
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The first two terms correspond to the present experimental constraints from solar neutrino 
oscillation and KamLAND summarized in eq.(6) *. In the last term, we assume that the planned 
future reactor experiments [25] should measure sin^ 2^rct with the expected uncertainty of 0.01. 



4 Octant degeneracy and the T2KK experiment 

In this section, we show the potential of the T2KK experiment for solving the octant degeneracy 
and investigate the role of the far detector and the future reactor experiments. We show in 
Fig.l the minimum Ax^ as a function of sin^ ^^tm expected at the T2KK experiment after 5 
years of data taking (5 x 10^^ POT). The event numbers are calculated for a combination of 
3.0° OAB at SK and 0.5° OAB at L = 1000km for the following parameters : 

{ml - m^)"'P"* = 2.5 x lO^^eV^ (normal hierarchy) , (27a) 

{ml - mlY'^P'"' = 8.2 x lO^^eV^ , (27b) 

sin2 2^^^P,"* = 0.10, (27c) 

sin' 2^X* = 0-83 , (27d) 

4T* = 0°,±90°,180°, (27e) 

^input ^ 3.0g/cm3 for L = lOOOkm , (27f ) 

^input ^ 2.8g/cm3 for SK . (27g) 

In the left-hand figure of Fig.l, we show the cases for the input values sin' =0.35 (a), 

0. 40 (b), 0.45 (c) and in the right-hand figure, for sin' ^i°P^* =0.55 (d), 0.60 (e), 0.65 (f). In 
each cases, the fit has been performed by surveying the whole parameter space. We find from 
Fig.l, that the octant degeneracy can be solved by T2KK experiment when sin' 2^atm = 0.91 

1. e., between sIu'^atm = 0.35 and 0.65 at 4(7. For sin'2^^ATM = 0.96 the degeneracy between 
sin'^ATM = 0.40 and 0.60 can be resolved with Ax' > 7, or at 2.6cr. However, it is difficult to 
solve the octant degeneracy for sin' 2^atm = 0.99, between sin' ^atm = 0.45 and 0.55. 

In the left-hand figure of Fig.l, the minimum Ax' for 6 = 90° is larger than those for the 
other CP phases. In the right-hand figure, the minimum Ax' is also largest aX 6 = 90°. There, 
however, the minimum Ax' for S = —90° is slightly larger than those for S = 0°, 180°. 

In order to explore the 5mns dependence of the capabihty of the T2KK experiment to solve 
the octant degeneracy, we show in Fig. 2 contours of the minimum Ax' in the whole space 
of sin'2^j,X* and S'^p^K The event numbers are calculated for various sin'2^j^P^* and S'^^f 
values in each figure, with sin' 6'^,^* = 0.40 for (a) and (c), or sin' 6^^™* — 0.60 for (b) and 
(d). The other model parameters are set as in eq.(27). Figs. 2 (a) and (b) are for the normal 
hierarchy, m\ — m{ — 2.5 x 10~^eV', and Figs. 2 (c) and (d) are for the inverted hierarchy, 
ml — ml — —2.5 x 10~^eV'. In performing the fit, all the 16 parameters (6 model parameters 
and 10 normahzation factors) are varied freely under the following constraints: sin' 6%^^ > 0.5 



*The most recent results, eq.(6), are slightly different from our inputs. Because our analysis is not sensitive 
to the difference, we use these values for the sake of keeping the consistency with our previous studies [8, 9]. 




Figure 1: Minimum Ax^ of the T2KK experiment as a function of sin^ ^atm- The event numbers 
are calculated for a combination of 3.0° OAB at SK and 0.5° OAB at L = 1000km with 100 
kton water Cerenkov detector, after 5 years running (5 x 10^^ POT). The input parameters are 
chosen as in eq.(27). In the left-hand figure, sin^ 6^^^* = 0.35 (a), 0.40 (b), 0.45 (c) and in the 
right-hand figure, sin^ ^ Q (^^-^^ q qq (^^^^ q (^^y 

for (a) and (c), sin^ 6'^!^j^ < 0.5 for (b) and (d), (m^-m^)*^* > for (a) and (b), (m^-m^)*^* < 
for (c) and (d). From Figs. 2 (a) and 2(c), we find that sin^ ^'^tm — O-^O can be distinguished from 
sin^ 9f^^ > 0.5 at Ax^ > 9 (4) for sin^ 26'™p^* ^0.12 (0.09) when the normal (inverted) hierarchy- 
is realized. Figs. 2(b) and 2(d) show that the octant degeneracy can be solved at Ax^ > 9 for 
sin^ 26'i,"Pr"* ^0.12 (0.14) when sin^ 6'^^^^ = 0.60 for the normal (inverted) hierarchy. 

It is found in Figs. 2 (a) and 2(b) that the minimum Ax^ is highest around = 90° 

confirming the trend observed in Fig.l. We find from Figs. 2(c) and 2(d) that the same trend 
holds even when the neutrino mass hierarchy is inverted. In all the four plots of Fig.2, we 
recognize a high plateau around ^^ns* — 90° ^ lower plateau around 5J^^g* = —90°. We 
can understand the trend by using the approximate expression of the z/^ Uf, transition 
probability, eq.(13). We first note that the Vf. oscillation probability is proportional to 

(1 + g)(l + A^^ sin^ 26'hct around the oscillation maxima, |Ai3 = {2n + l)vr|. When q = —0.2, 
the z/g appearance rate is proportional to 0.8(1 + A'^). In order to reproduce the same rate for 
q = 0.2, we should find a parameter set that makes the factor 1 + 40% smaller than its input 
value, up to the uncertainty in sin^26'j^cT5 which is assumed to be 0.01/sin^ 26^^^* in eq.(26). 
This cannot be achieved for 5j^f^g* = 90°, because the input value of 1 + takes its minimum 
value. On the other hand, if ^mns* — ~90° the input value of 1 + A"^ is large and it can be 
reduced siginificantly by choosing 6^-^^ = 90° in the fit. This explains why the minimum Ax^ 
is larger around ^^Ng* = 90° than that around ^i^Jg* = —90° when sin^ ^^tm = 0.4 in Figs. 2(a) 
and 2(c). When q = 0.2, the same argument tells that we cannot compensate for the large 
input value of (1 + g)(l + A'^) for 6^^^^^ = —90°. This indeed explains the lower plateau around 
^MNs* = ~90° observed in 2(b) and 2(d). The cause of the higher plateau around 5mns* = 90° 
in these figures for sin^ 6'atm = 0.6 is more subtle. When 5™^"* = 90°, (1 + ^^ynput ^^i^gg j^g 
smallest value, and the reduction in (1 + q) from (1 + g)™?"* = 1 + 0.2 to (1 + q)^^ = 1 — 0.2 
can be compensated for by making (1 + A^)^^ larger by choosing 5^*j^g ~ —90°. This, however. 
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Figure 2: The potential of the T2KK experiment to solve the octant degeneracy with the same 
OAB combination of the Fig.l. In each figures, the event numbers are obtained for the model 
parameters at various sin^ 2^^^"*. K^T and sin^ ^jX* = 0.40 for (a,c), or sin^ ^jX* = 0.60 
for (b,d), with the normal (a,b) or the inverted (c,d) hierarchy. The other input parameters 
are same as those of Fig.l. All the parameters are taken freely in the fit under the constraint 
sin^ 6'^!^^ > 0.50 (a,c), or sin^6'^!^j^ < 0.50 (b,d). The resulting values of minimum Ax^ are 
shown as contours for Ax^ =4, 9, 16. 

necessarily makes the coefficient of lAisl/vr in eq.(13b) have the wrong sign, and hence the ratio 
of the first peak (|Ai3| ~ tt) and the second peak (|Ai3| ~ Svr) cannot be reproduced. The 
higher plateau around 5™^"* = 90° in Figs.2(b) and 2(d) for sin^ = 0-6, and the lower 

plateau around = -90° in Figs.2(a) and 2(c) for sin^ ^jX* = 0.4 can be explained as 

above. 

We show in Fig. 3 the allowed region of sin^ ^^tm and sin^ 26'kct by the T2KK experiment. 
The event numbers are generated at = 0° and sin^ 2^^=^p,"* = 0.10 for sin^ ^^X* = 0-35, 

0.40, 0.45, 0.50, 0.55, 0.60, 0.65 from the 1st to the 7th row. The other input parameters are 
the same as in eq.(27). The allowed regions in the plane of sin^ 6'atm and sin^ 2^kct are shown 
by the Ax^ = 1, 4, 9 contours depicted as solid, dashed, and dotted lines, respectively. In 
the left-hand-side plots, (a), the constraint on sin^ 26'j,ct from the future reactor experiment is 
kept in the Ax^ function. On the other hand, in the right-hand-side plots, (b), the external 
constrant on sin^ 2^hct is removed from the Ax^ function in eq.(26). Comparing Figs. 3(a) and 
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Figure 3: The capability of the T2KK experiment for constraining the sin^ 6'atm and sin^ 2^^rct- 
Allowed regions in the plane of sin^ 6'atm and sin^ 29^ct are shown for the same T2KK set up 
Fig.l. In each figure, the event numbers are generated at 5™Pg* = 0° and sin^ 26^^"* = 0.10, for 
sin^^i'^P^* = 0.35, 0.40, 0.45, 0.50, 0.55, 0.60, 0.65 from the 1st to the 7th row. The other input 
parameters are the same as in Fig.l. In the left-hand-side plots, (a), we keep the constraint 
on sin^ 2^RCT form the future reactor experiment, whereas in the right-hand-side plots, (b), we 
remove the external constraint on sin^ 2^rct in eq.(26). The Ax^ = 1, 4, 9 contours are shown 
by the solid, dashed, and dotted lines, respectively. 
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Figure 4: The same as Fig.2, but with 5 times larger exposure (25 x 10^^ POT). 
3(b), we find that the mirror solution around 



sin^ 26 



fit 



sin^ 2fe"* 
1-q ''''^ 



(2J 



cannot be excluded without the information from the future reactor experiment. 

Before closing the section, we examine the impact of upgrading the J-PARC beam intensity 
be a factor of 5 [26] on the resolution of the octant degeneracy. Such an upgrade is desirable 
especially if the neutrino mass hierarchy is inverted, because the octant degeneracy between 
sin^ 6'atm = 0.4 and 0.6 cannot be resolved at 3a unless 6mns — 90°; see Figs. 2(c) and 2(d). 

We show in Fig. 4 the same contour plots as in Fig.2, but with 5 times larger exposure 
(25 X 10^^ POT). It is found that the degeneracy between sin^ 6'atm = 0.4 and 0.6 can now 
be resolved at 3cr level for sin^ 26'rct > 0.08(0.09), when the hierarchy is normal (inverted). 
Comparing Fig.2 and Fig. 4, however, we find that the sensitivity does not improved as much 
as we would hope with 5 times higher statistics. The minimum Ax^ value does not grow by 
a factor 5, because the capability of resolving the octant degeneracy is now dictated by the 
accuracy of the external constraint on sin^ 26'rct from the future reactor experiment. 



sin^ 2e„c 



0.01 



sin^ 26'rct ' 



(29) 



which we assume in eq.(26). The fractional uncertainty of sin^ 2^^rct is 10% for sin^ 29^ct = 0.1, 




Figure 5: Minimum Ax^ as a function of sin^ ^Ttm ^'^^ ^^e same T2KK setting as in Fig.l. In 
each figure, the event numbers are calculated for the parameters of eq.(27) with 5™^"* = 0° 
(a), 90° (b), 180° (c), —90° (d), under the normal hierarchy, and the fit has been performed 
by assuming the inverted hierarchy. The solid line gives the minimum Ax^. The open circle 
(square) denotes the minimum value of Ax^ when the sign of q^^P'^^q^^ is positive (negative). 

but it is 17% for sin^ 2^rct = 0.06. If sin^ 29^ct turns out to be even smaller, the fractional 
error grows and the mirror solution eq.(28) can no more be resolved. If sin^ 2^j^ct turns out 
to be smaller than 0.06, further reduction of its error in the future experiments with reactor 
and/or the beta beam [27]. 



5 Mass hierarchy and the octant degeneracy 

In this section, we examine the effect of the octant degeneracy on the capability of the T2KK 
experiment to determine the neutrino mass hierarchy pattern. 

Figures shows the minimum Ax^ as a function of sin^ ^Ttm ^^e T2KK experiment to 
determine the mass hierarchy pattern with the same OAB combination of Fig.l. In each figure, 
the event numbers are calculated for 5jjPf = 0° (a), 90° (b), 180° (c), and -90° (d), when 
the normal hierarchy is realized. The other parameters are listed in eq.(27). The fit has been 
performed by surveying the whole parameter space by assuming the wrong hierarchy. The solid 
line gives the minimum Ax^. The open circle denotes the minimum Ax^ when the sign of 
g,mput^fit^ or that of (1 — 2 sin^ ^Ttm){^ ~ 2 sin^ ^^m)^ is positive, whereas the open square gives 



the minimum A^^ when q^^p^^q^^ is negative. 

When g^* takes the same sign as g^^^P"*, sin^^f!^^ ~ sin^ ^JJ^^* is favored, and the reduction 
of the — > oscillation amplitude (1 + A^) in eq.(13b) for the inverted hierarchy, A13 ~ 
—71, cannot be compensated for in the two detector experiment [8, 9]. Because the Ue 
rate is proportional to sin^^^,^*, the resulting increase in the discrepancy leads to the linear 
dependence of the minimum Ax^ on sin^ ^"tm* observed for the open circle points. On the 
other hand, when sin^^J^pJJ* < 0.5 (g™P"t < 0), it is possible to compensate for the reduction 
of the 1 + factor of the Ug oscillation amplitude by choosing q^^ ~ _ginput^ since 

sin^ O^^^i = sin^ ( 1 + + g"^P'it) > 0.5. This explains why the open square points 

for g"ipitgfit < gives the lowest Ax^ for sin^ GTtm < When sin^ ^atm* is significantly 

lower than 0.5, however, the enlargement factor of sin^ ^atm/ ^^^'^ ^atm = (1 + Q'^*)/ (1 + 
overshoots the reduction due to the matter effect, especially for the z/^ rate at SK when the 
matter effect is small. When ^^^s* — ~90°, shown in Fig. 5(d), this reduction in the minimum 
Ax^ by using the octant degeneracy is most significant because the overshooting of the — > Ve 
rate can be partially compensated by choosing sin(5^*j^g > 0. 

In Fig. 6, we show the minimum Ax^ as a function of sin^ ^^m* when the neutrino mass 
hierarchy is inverted. The gradual increase of the Ax as sin^ ^atm* grows can also be seen for the 
open circle points where the fit is restricted to the parameter space that satisfies gStg^P"* > 0. 
This refiects the increase of the v^. event rate as sin^ ^^m* increases, independent of the 

hierarchy pattern. On the other hand, the minimum Ax^ for the parameter space of (f^q^^^^^ < 
0, plotted by open squares, gives the lowest Ax^ when sin^ ^^m* > O-^- This is because 
the reduction of the v,j, rate in the fit, which is proportional to sin^ ^a™/ sin^ ^^m* = 

(1 + g^*)/(l + g™P"*) < 1, for g*"* < < g'"?"*, can compensate for the reduction due to the 
matter effect when the hierarchy is inverted. The reduction of the minimum Ax^ due to the 
octant degeneracy is strong at sin^ ^^m* > 0.5 for the inverted hierarchy, and the increased 
sensitivity to the mass hierarchy pattern for large sin^ ^atm is lost when it is inverted. 

In Fig. 7, we show the capabihty of the T2KK experiment to determine the mass hierarchy 
pattern as contour plots of the minimum Ax^ value on the parameter space of sin^ 2^™p"* 
and (^MNs*- -'-^ GSiCh. figure, the input date are calculated for the model parameters at various 
sin^2e'^^^^ and values, with sin^ = 0.40 in (al, bl), sin^ ^jX* = 0.50 in (a2, b2), or 
sin^ ^S?M* — 0-60 in (a3, b3). The left-hand figures (al, a2, a3) are for the the normal hierarchy, 
and the right-hand figures (bl, b2, b3) arc for the the inverted hierarchy. The other input 
parameters are the same as those of Fig. 5, and in eq.(27). All the fit parameters are varied 
freely to minimize the Ax^ function, under the constraint of the opposite mass hierarchy. The 
resulting values of minimum Ax^ are shown as contours for Ax^ = 9, 16, 25. The contours of 
Figs.7(a2) and 7(b2) are identical to those of Fig. 6 of Ref. [9], which we copy for the purpose 
of comparison. 

It is clearly seen in Fig. 7 that the main feature of the T2KK ability for the mass hierarchy 
determination at sin^ 2^^^,^* = 0.96, sin^ 6'atm = 0.4 (al, bl) or 0.6 (a3, b3), arc not much 
different from those at sin^ 2^^^* = 1.0 (a2, b2), such as the fact the minimum Ax^ around 
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Figure 6: The same as Fig. 5, but when the inverted hierarchy is reahzed in nature and the fit 
is performed by assuming the normal hierarchy. 

(5™Pg' ~ 0° is smaller than that around 5™^g* ~ 180°. Close examination of Fig. 7, however, re- 
veals the foUowings. In case of the normal hierarchy, m\—m\ > 0, the minimum Ax^ grows with 
growing sin^ ^^m*; the mass hierarchy can be determined at 3a level for sin^ 26'™^"* ^ 0.07 
if sin^^jfP^t = 0.40 (al), whereas the same holds for sin^ 2^j,'g'^"* ^0.04 if sin^^^P^* = 0.60 (a3). 
This is mainly because the z/^ Ve event rate grows with sin^ G^atm ^"^^ because the presence 
of the octant degeneracy does not disturb the measurement significantly, as can be seen from 
the open circle points in Fig.5. In contract, no significant improvement in the hierarchy dis- 
crimination power is found for sin^ (^""atm > 0-^ in case of the inverted hierarchy. This is because 
the octant degeneracy between sin ef^^ = 0.6 and sin^ = 0.4 allows us to compensate for 
the matter effect reduction of the z/^j Ve rate. We find that the best hierarchy discrimination 
is achieved at sin^6'^j^* = 0.5 for all 5mns values, confirming the trends of Fig. 6. 

Summing up, the T2KK two detector experiments can resolve the mass hierarchy pattern in 
the presence of the octant degeneracy. If the hierarchy is normal, the discriminating power grows 
with increasing sin^ (^^atm- ^^e other hand, if the hierarchy is inverted, the discriminating 
power reduces both at sin 

^iX* < 0.5 and at sin^^jfP^f > 0.5: it reduces at sin ^ATM < 0.5 
because of the lower rate of the — > events, while it reduces at sin^ 6'atm > 0.5 because of 
the octant degeneracy. 
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Figure 7: The potential of the T2KK experiment to determine the neutrino mass hierarchy. 
The input data calculated for the normal hierarchy (al, a2, a3) or for the inverted hierarchy 
(bl, b2, b3), and the fit has been performed by assuming the wrong hierarchy. In each figure, 
the minimum Ax^ is obtained for the input data calculated at various (sin^ 26'™^"*, 5mns*) with 
sin^^i^P^t = 0.40 for (al, bl), sin^ 0^p^* = 0.50 for (a2, b2), or sin^ 0jf p^* = 0.60 for (a3, b3). 
The other input parameters are the same as those of Fig.l. The contours of the minimum Ax^ 
are shown for Ax^ = 9, 16, 25. 



180 



- 



-90 

































- 
- 




























- 

- 


- 




























- 


- 




























■ 


















































% 

































































































































































































-180 
0.00 



0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 
sin^2e„cT 



































(c) 




















































































































































































































































































































1 





















0.02 0.04 



0.06 0.08 

sin=^2e„cT 



0.10 0.12 0.14 



-90 



180 

































\") 






































< 


1 






















































































































































































































































































DO 


0.02 


0.04 


0.06 0.08 
sin220RCT 


0.10 


0.12 


0.14 



-90 



:(d): 



u'0 



0.00 0.02 0.04 



0.06 0.08 
sin22e„cT 



0.10 0.12 



0.14 



Figure 8: The allowed region in the plane of ^Yc?1d^c^ and 5mns by the T2KK set up of Fig.l, 
when sin^ ^atm* — 0.40 and the hierarchy is normal. The input values of sin'^ '2B^q^ and 5mns; 
sin^ 2^j,°Pr"* = 0.06 or 0.10, b'Zf = 0° (a), 90° (b), 180° (c), and -90° (d), are denoted by solid 
blobs in each figure and the other input parameters are listed in eq.(27). The la, 2a, and 3a 
contours are shown by solid, dashed, and dotted lines, respectively. The thick (thin) lines are 
for sin^ 2^^^^"* = 0.10 (0.06). There is no allowed region within 3a when the inverted hierarchy 
is assumed with fit. 



6 CP phase and the octant degeneracy 

In this section, we investigate the relation between the CP phase measurement and the octant 
degeneracy. 

Figure8 shows the potential of the T2KK experiment for measuring sin^ 20rct and 5mns when 
sin^ ^^M* = 0.40 and the hierarchy is normal, rrig — > 0. The input values of sin^ ^O^ct and 
5mns are denoted by solid blobs in each figure, and the la, 2a, and 3a allowed regions are shown 
by solid, dashed, and dotted lines, respectively. The thick contours are for sin^ 2^^^* = 0.10, 
and the thin contours are for sin^ 2^^X"* = 0-06. = 0° in Fig.8(a), 90° in (b), 180° in (c), 

and —90° in (d). There is no additional allowed region within 3a when the inverted hierarchy 
is assumed in the fit, in accordance with Fig.7(al). 

When we compare the contours of Fig. 8 with the corresponding ones in Fig. 8 of Ref. [9] 
for sin^ (^^atm = 0-5, we can clearly identify the islands due to the octant degeneracy. When 
sin^ 2^™P"* = 0.10, the thick contours shows no island at 5mns* — 90°' but 3a islands appear at all 
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Figure 9: The same as Fig. 8, but for sin^ ^Ttm = 0.60. 



the other S^^^^*' cases, which is consistent with the result of Fig. 2(a). In case of sin^ (^rct^ = 0.06, 
the contours have 3a islands for all ^mns* cases and a clear 2a island at S^^^^^ = —90°, again 
consistent with Fig.2(a). Close examination of the location of the islands reveals that their 
center is at around sin^ 26'|J;^ = sin^ 26'™(?"* (0.4)/(0.6), as expected by eq.(28). The existence 
of the islands due to the octant degeneracy hence reduces our capability of measuring sin^ 2^rct 
significantly. 

It is remarkable that the octant degeneracy does not jeopardize the T2KK capability of 
determining the CP phase, 6mns- the islands in Fig. 8 have the 5mns values consistent with its 
input values. This is because the coefficients of both sin 5iy[i,js and cos (^m^s 

in the — > 

oscillation probability is not sensitive to the octant degeneracy, as explained in eq.(14). The 
results we found in Fig. 8 confirms the validity of our approximation for the T2KK experiments. 

Figure9 also shows the potential of the T2KK experiment for measuring sin^ 26'rct and 6mns, 
but for sin^ ^Ttm — 0.60. There is no allowed region within 3a when the inverted hierarchy is 
assumed in the fit, as can be seen from Fig.7(a3). 

Comparing the contours of Fig. 9 with the corresponding ones in Fig.8 of Ref. [9] for sin^ (^Ttm — 
0.5, we can again identify the islands due to the octant degeneracy. When sin^ 26*^"* = 0.10, 
the thick contours shows no island at (JJ^^g* = ±90°, but 3a islands appear at the other ^^ns* 
cases. In case of sin^ ^^bct^ — 0.06, the contours have 3a islands for all ^^^s* cases, but there 
is no 2a island. These results are consistent with the result of Fig.2(a). The location of the 
center of the islands is at around sin^ 26'||;^ = sin^ 26'^'^p"* (0.6)/(0.4), as expected by eq.(28). 



The capability of measuring sin^ 2^rct for sin^ ^JSm* = 0.60 is also reduced by the octant de- 
generacy, which is the same as that for sin^ ^Ttm — 0-40. However, the octant degeneracy does 
not disturb the T2KK capability of determining the CP phase. 

7 Summary 

There are three types of ambiguity in the neutrino parameter space. The first one is from 
the sign of the larger mass-squared difference, which is related to the mass hierarchy pattern. 
The second one is from the combination of the unmeasured parameters, the leptonic CP phase 
(5mns) and the mixing angle ^rct, which resides at the upper-right corner of the MNS matrix 
[1]. In the previous studies [8, 9], we showed that the idea [10] of placing a lOOkton-level 
water Cercnkov detector in Korea along the T2K neutrino beam at L = 1000km, the T2KK 
experiment [13] can solve these ambiguities. The last ambiguity is in the value of the ^atm, which 
dictates the atmospheric neutrino observation [2] and the long base-line neutrino oscillation 
experiment [3, 4]. If the mixing angle ^^atm, is not 45°, there is a two fold ambiguity between 
"^atm" and "90° - 6 

atm"; the octant degeneracy [16]. 

In this paper, we focus on the physics potential of the T2KK experiment for solving the 
octant degeneracy. In our scmi-quantitatively analysis, wc follow the strategy of Rcf. [8, 9] 
where wc adopt SK as a near side detector and postulate a 100 kton water Ccrenkov detector 
at L = 1000km, and the J-PARC neutrino beam orientation is adjusted to 3.0° at SK and 0.5° 
at the Korean detector site. 

If the value of sin^ 2^atm is 0.99, which is 1% smaller than the maximal mixing, the value 
of sin ^atm is sin ^atm — 0.45 or 0.55, which differ by 20%. Therefore, we also investigate the 
impacts of the octant degeneracy on the physics potential for the mass hierarchy determination 
and the CP phase measurement by T2KK, because the leading term of — > Ue oscillation 
probability is proportional to sin^ ^atm, not sin^ 2^atm- 

When we include the constraint for the value of sin^^Rcx, which will be obtained from the 
future reactor experiments [25]. the octant degeneracy between sIu^^atm = 0.40 and 0.60 can 
be resolved at 3a level for sin^ 2^j^ct > 0.12 (0.08) after 5 years exposure with 1.0 (5.0) x 10^^ 
POT/year, if the hierarchy is normal, see Fig.l and Fig. 2. We find that the contribution from 
the second maximum of the — > oscillation probability at the far detector (L = 1000km) 
plays an important role for solving the octant degeneracy. It is also found that the octant 
degeneracy cannot be solved without the contribution from future reactor experiments. 

We also investigate the impact of the octant degeneracy in the determination of the mass 
hierarchy pattern. The T2KK power of resolving the mass hierarchy pattern is proportional 
to the value of sin^ (^^tm fo^ the normal hierarchy, see Fig. 5, because the z/^ rate is 

proportional to sin^ ^^J^j^'J' . When the mass hierarchy is normal, we can determine the mass 
hierarchy at 3(j level for ^\Y?2Q^^lf- ^0.07 if sin^ 6'j3Pf = 0.40, Fig.7(al), whereas the same 
holds for ^m^^B'^lf' ^0.04 if sin^ ^i'^P.f = 0.60, Fig.7(a3). On the other hand, if the hierarchy 
is inverted, sin^ ^atm — 0.5 is found to be the optimal case for the mass hierarchy determination. 



see Fig.6, because of the lower rate of the — > events for sin^ ^^^m* < 0.5 and the octant 
degeneracy for sin^ ^atm > 

Finally, we check the effect of the octant degeneracy for the CP phase measurement, see 
Fig. 8 and Fig.9. The CP phase can be constrained to ±30° at la level for sin^2^ATM = 0.96, 
even if we cannot distinguish between sin ^^ATM — 0.4 or 0.6. The error does not increase from 
that for sin^ ^atm = 0.5, because the coefficients of both sine and cosine term of Smns in the 
oscillation probability are not sensitive to the octant degeneracy. 
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Note Added 

When finalizing the manuscript, we learned that a similar study has been performed by T. Ka- 
jita, et al. [28]. 
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